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Abstract 

We study the production of the Higgs bosons predicted in the Minimal Su- 

persymmetric extension of the Standard Model (h° , H° , A , H^), with the 

reactions e + e~ — ► bbh (H° , A ), and e + e~ — ► r~ 9 T H + ,t + u t H~ , using the 

helicity formalism. We evaluate cross section of h°,H°,A° and in the 

limit when tan f3 is large. The numerical computation is done considering 

two stages of a possible Next Linear e + e~ Collider: the first with ^/s = 500 

GeV and design luminosity 50 fb" 1 , and the second with y/s = 1 TeV and 

luminosity 100-200 fb' 1 . 
PACS: 14.80.Cp, 12.60.Jv 
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I. INTRODUCTION 



Higgs bosons [|I| play an important role in the Standard Model (SM) 0; they are respon- 
sible for generating the masses of all the elementary particles (leptons, quarks, and gauge 
bosons). However, the Higgs-boson sector is the least tested one in the SM. If Higgs bosons 
are responsible for breaking the symmetry from SU{2) L x U(l)y to U(1)em, it is natural 
to expect that other Higgs bosons are also involved in breaking other symmetries. One of 
the more attractive extensions of the SM is Supersymmetry (SUSY) ||, mainly because of 
its capacity to solve the naturalness and hierarchy problems while maintaining the Higgs 
bosons elementary. 

The theoretical frame work of this paper is the Minimal Supersymmetric extension of 
the Standard Model (MSSM), which doubles the spectrum of particles of the SM and the 
new free parameters obey simple relations. The scalar sector of the MSSM requires two 
Higgs doublets, thus the remaining scalar spectrum contains the following physical states: 
two CP-even Higgs scalar (h° and H°) with m h o < m H o, one CP-odd Higgs scalar (A ) and 
a charged Higgs pair (H^, whose detection would be a clear signal of new physics. The 
Higgs sector is specified at tree-level by fixing two parameters, which can be chosen as the 
mass of the pseudoscalar m^o and the ratio of vacuum expectation values of the two doublets 
tan (3 = Vijv i , then the mass m^o , m#o and m H ± and the mixing angle of the neutral Higgs 
sector a can be fixed. However, since radiative corrections produce substantial effects on 
the predictions of the model ||, it is necessary to specify also the squark masses, which are 
assumed to be degenerated. In this paper, we focus on the phenomenology of the neutral 
CP-even and CP-odd scalar (h°,H°,A°) and charged {H^). 

The search for these scalars has begun at LEP, and current low energy bound on their 
masses gives m h o, m A o > 90 GeV and m H ± > 120 GeV for tan /3 > 1 @. At e + e~ colliders the 
signals for Higgs bosons are relatively clean and the opportunities for discovery and detailed 
study will be excellent. The most important processes for the production and detection of 
the neutral and charged Higgs bosons h°, H°, A and H^, are: e + e~ — > Z* — > h°, H° + Z°, 
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e+e" -> Z* -> + e+e" -> -> vv + h°, H° (the later is conventionally 

referred to as WW fusions), and e + e~ — > H + H~ 0; precise cross-section formulas appear 
in Ref. ||. The main decay modes of the neutral Higgs particles are in general bb decays 
(~ 90%) and t + t _ decays (~ 10%) which are easy to detect experimentally at e + e~ colliders 
PHTT|. Charged Higgs particles decay predominantly into tv t and tb pairs. 

In previous studies, the two-body processes of the neutral and charged Higgs bosons 
e+e" -> h°(H°) + Z°, e+e" -> h°(H°) + A and e+e" -> H + H~ have been evaluated § 
extensively. However, the inclusion of three-body process with heavy fermions /, e + e _ — > 

and e + e _ — ► t~v t H + ,t + v t H~ at future e + e~ colliders 

energies is necessary in order to know its impact on the two-body mode processes 

and also to search for new relations that could have a cleaner signature of the Higgs bosons 
production. 

In the present paper we study the production of SUSY Higgs bosons at e + e~ colliders. 
We are interested in finding regions that could allow the detection of the SUSY Higgs bosons 
for the set parameter space (m^o — tan/3). We shall discuss the neutral and charged Higgs 
bosons production bbh°(H°, A ), and t~v t H + , t + v t H~ in the energy range of a future e + e~ 
collider ]T^-0 for large values of the parameter tan/5, where one expects to have a high 
production. Since the coupling h°bb is proportional to sin a/ cos/5, the cross-section will 
receive a large enhancement factor when tan/3 is large. Similar situation occurs for H°, 
whose coupling with bb is proportional to cos a/ cos (3. The couplings of A with bb and of 
with t u T , t + u t are directly proportional to tan f3, thus the amplitudes will always grow 
with tan (3. We consider the complete set of Feynman diagrams at tree-level and use the 
helicity formalism |19^p5[ for the evaluation of the amplitudes. Succinctly, our aim in this 
work is to analyze how much the results of the Bjorken Mechanism [Fig. 1, (1.4)] would 
be enhanced by the contribution from the diagrams depicted in Figs. 1.1-1.3, 1.5 and 1.6 
in which the SUSY Higgs bosons are radiated by a b(b) quark. For the case of the charged 
Higgs bosons the two-body mode [Figs. 3.1 and 3.4] would be enhanced by the contribution 
from the diagrams depicted in Figs. 3.2, 3.3, and 3.5, in which the charged Higgs boson is 
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(ffh°JfH°JfA% [12,13 



radiated by a t~u t {j + u r ) lepton. 

Recently, it has been shown that for large values of tan f3 the detection of SUSY Higgs 
bosons is possible at FNAL and LHC f2H[| . In the papers cited in Ref. |26[ the authors 



calculated the corresponding three-body diagrams for hadron collisions. They pointed out 
the importance of a large bottom Yukawa coupling at hadron colliders and showed that the 
Tevatron collider may be a good place for detecting SUSY Higgs bosons. In the case of the 
hadron colliders the three-body diagrams come from gluon fusion and this fact makes the 
contribution from these diagrams more important, due to the gluon abundance inside the 
hadrons. The advantage for the case of e + e~ colliders is that the signals of the processes are 
cleaner. 

This paper is organized as follows. We present in Sec. II the relevant details of the 
calculations. In Sec. Ill we evaluate cross section for the processes e + e~ — > bbh°(H°,A ) 
and e + e~ — > t~u t H + , t + v t H~ at future e + e~ colliders. Finally, Sec. IV contains our 
conclusions. 

II. HELICITY AMPLITUDE FOR HIGGS BOSONS PRODUCTION 

When the number of Feynman diagrams is increased, the calculation of the amplitude is 



a rather unpleasant task. Some algebraic forms |27j can be used in it to avoid manual calcu- 
lation, but sometimes the lengthy printed output from the computer is overwhelming, and 
one can hardly find the required results from it. The CALKUL collaboration p8| suggested 
the Helicity Amplitude Method (HAM) which can simplify the calculation remarkably and 
hence make the manual calculation realistic. 

In this section we describe in brief the evaluation of the amplitudes at tree-level, for 
e + e- -> bbh°(H°, A ) and e + e~ -> r-jy T H + , r + u T H- using the HAM This method is 

a powerful technique for computing helicity amplitudes for multiparticle processes involving 
massless spin- 1/2 and spin-1 particles. Generalization of this method that incorporates 



massive spin- 1/2 and spin-1 particles, is given in Ref. K5j. This algebra is easy to program 
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and more efficient than computing the Dirac algebra. 

A Higgs boson h°, H°, A , and H ± can be produced in scattering e + e~ via the following 
processes: 

e+e~ -> bbh°, (1) 

e+e" -> (2) 

e+e" -> (3) 

e+e" -> t-u t H + , t + v t H-. (4) 

The diagrams of Feynman, which contribute at tree-level to the different reaction mech- 
anisms, are depicted in Figs. 1-3. 

However those diagrams with exchange of Higgs bosons instead of gauge bosons (photon 
or Z°) have been neglected because of the smallness of the Higgs fermion coupling. Using 
the Feynman rules given by the Minimal Supersymmetric Standard Model (MSSM), as 
summarized in Ref. ||, we can write the amplitudes for these reactions. For the evaluation 
of the amplitudes we have used the spinor-helicity technique of Xu, Zhang and Chang |2H 



(denoted henceforth by XZC), which is a modification of the technique developed by the 
CALKUL collaboration |28|| . Following XZC, we introduce a very useful notation for the 
calculation of the processes (l)-(4). The complete formulas of the processes (l)-(3) are given 
in Ref. [ 13| , while for the case of process (4), are given in Ref. ||14j| . We are going to use the 



same notation and the formulas given in Hi3| , |i4|1 and do not reproduce them here. 

After writing down the Feynman diagrams corresponding to a given amplitude Ai an 
usually proceeds to derive an analytic expression for the cross section Sl-M) 2 ; with an 
appropriate spin and/or color sum or average. The result, which is usually a function of 
Minkowsky products of the particle four-momenta, is then evaluated numerically at phase- 
space points in the region of interest. 

It is clear that any method which evaluates A4 instead of J2 \-M\ 2 will eventually become 
superior to the standard approach. Indeed, several authors have stressed this point and 
proposed alternatives. 



We want to argue that it is useful to employ in the evaluation of the amplitude not vector 
products like p\ ■ p 2 but rather spinor products like u(pi)u(p 2 ). 

After the evaluation of the amplitudes of the corresponding diagrams, we obtain the cross- 
sections of the analyzed processes for each point of the phase space using Eqs. (12)-(17), 
(23)-(28) and (11)-(15) [0,0 by a computer program, which makes use of the subroutine 
RAMBO (Random Momenta Beautifully Organized) |29j . The advantages of this procedure 
in comparison to the traditional "trace technique" are discussed in Refs. Hl9|-[25[| . 

We use the Breit-Wigner propagators for the Z°, h°, H°, A and H ± bosons. For the 
SM parameters we adopted the following: m;, = 4.25 GeV, rrit = 175 GeV, m T = 1.78 
GeV, m v = 0, m z o = 91.2 GeV, T z o = 2.4974 GeV, sin 2 6 W = 0.232, which are taken as 
inputs. The widths of h°, H°, A and are calculated from the formulas given in Ref. ||. 
In the next sections we present the numerical computation of the processes e + e~ — > bbh, 
h = h°, H°, A and e+e~ -> t u t H^, t+v t H-. 



III. PRODUCTION OF THE MSSM HIGGS BOSONS AT NEXT GENERATION 

POSITRON-ELECTRON COLLIDERS 

In this paper, we evaluate total cross section of neutral and charged MSSM Higgs bosons 
at next generation e + e~ colliders, including three-body mode diagrams [Figs. 1.1-1.3, 1.5, 
and 1.6; Figs. 2.1-2.3, 2.5 and 2.6; Figs. 3.2, 3.3, and 3.5] besides the dominant mode 
diagram [Fig. 1.4; Fig. 2.4; Figs. 3.1, and 3.4] consider two stages of a possible Next Linear 
e + e~ Collider: the first with a/s = 500 GeV and design luminosity 50 fb , and the second 
with = 1 TeV and luminosity 100-200 fb" 1 . We consider the complete set of Feynman 
diagrams (Figs. 1-3) at tree-level and utilize the helicity formalism for the evaluation of 
their amplitudes. In the next subsections, we present our results for the case of the different 
Higgs bosons. 
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A. Cross section of h° 



In order to illustrate our results on the production of the h° Higgs boson, we present 
graphs of the cross section as functions of A Higgs boson mass m^o, assuming m t = 175 
GeV, M~ = 500 GeV and tan (3 > 1 for NLC. Our results are displayed in Fig. 4, for the 
process at three-body e + e~ — > bbh° and for the e + e~ — > (A , Z°) + h° dominant mode. 

The total cross section at colliders energies of 500 GeV and 4 different values of the 
fundamental supersymmetry parameter tan/3, 6, 10, 30, 50 is of the order of 0.1 pb. We 
note from this figure, that the effect of the reaction bbhP is light more important that 
(A , Z°) + h°, for the interval 80 GeV< m A o < 110 GeV and tan/3 = 30,50. Nevertheless, 
there are substantial portions in which the discovery of the h° is not possible using either 
bbh° or {A°,Z°) + h°. 

For the case of y/s = 1 TeV, the results of the total cross section of the h° are shown 
in Fig. 8. It is clear from this figure that the contribution of the process e + e~ — > bbh° 
becomes dominant in the interval 80 GeV < m A o < 120 GeV. However, they could provide 
important information on the Higgs bosons detection. For instance, we present in the next 
section tables that illustrate the events number for the process e + e~ -> bbh°, with m^o = 100 
GeV, v/s = (500, 1000) GeV and tan/3 = 10, 30. 

B. Cross section of H° 

To illustrate our results regarding the detection of the heavy Higgs bosons H°, we give 
the total cross section for both processes e + e" — > bbH° and e + e~ — > (A , Z°) + H° in Fig. 
5 for v/i = 500 GeV. 

The total cross section for this case is of the order of 0.01 pb. In this figure, we observed 
that the effect of incorporate bbH° in the detection of the Higgs boson H° is more important 
than the case of two-body mode (A , Z°) + H°, because bbH° is dominant in all the interval 
of 100 GeV< m A o < 400 GeV. 



7 



For the case of = 1 TeV, the results are show in Fig. 9. In this case the three-body 
mode bb~H° is light more important that (A , Z°) + H° for 100 GeV< m A o < 450 GeV . In 
the next section are present tables with the events number for e + e~ — > bbH°. 

C. Cross section of ^4° 

For the pseudoscalar A , it is interesting to consider the production mode in bbA°, since 
it can have large a cross-section due to the fact that the coupling of A with bb is directly 
proportional to tan/3, thus will always grow with it. In Fig. 6, we present the total cross 
section for the process of our interest bbA°, at \fs = 500 GeV and tan/3 = 6, 10, 30, 50. 

We note that for tan/3 = 6, 10, 30, 50, both the bbA° mode and the (h°, H°) + A mode 
have cross sections of the order of 0.1 pb, and for m A o > 175 GeV the process bb A cover a 
major region in the parameters space (m A o , tan ($) . 

On the other hand, if we focus the detection of the A at y/s = 1 TeV, the panorama for 
its detection is more extensive. The Fig. 10 shows the contours lines in the plane (a T — m A o), 
to the production cross section of bbA°. The contours for this cross section correspond to 
tan/3 = 6, 10, 30, 50. The total cross section is of the order of 0.01 pb and the three-body 
mode e + e~ — > bbA° is more important that two-body mode e + e~ — > (h°, H°) + A in all the 
parameters space. 

D. Cross section of H ± 

Our results for the H + H~ scalars are displayed in Fig. 7, 11 for the processes at three- 
body e + e~ — > t~u t H + ,t + u t H~ and for the e + e~ — > H + H~ dominant mode. 

The total cross section for this reaction with y/s = 500 GeV and 4 different values of 
tan/3, 6, 10, 30, 50 are show in Fig. 7. The cross section is of the order of a T ~ 0.1 pb for 
m H ± 100 GeV and tan/3 large, while for m H ± m 200 GeV is of a T m 0.01 pb. We can 
see from this figure that the effect of the reactions t~v t H + and t + v t H~ is slightly more 
important than H + H~ for most of the (m A o —tan /3) parameters space regions. It is precisely 

8 



in this curve where the contribution of the processes at three-body is notable. Nevertheless, 
there are substantial portions of parameters space in which the discovery of the H ± is not 
possible using either H + H~ or t~u t H + and t + v t H~ . 

For the case of y/s = 1 TeV, the results are show in Fig. 11. In both cases (t~v t H + , 
t + i* t H~ and H + H ) the curves with tan/5, 6, 10, 30, 50 given <j T ~ 0.1 pb for m H ± ps 175 
GeV. While for m H ± m 350 GeV the cross section drop until 0.001 pb^ 1 for tan/3 = 6. 
These cross sections are small, however, the contribution of the processes at three-body is 
slightly more important than the process at two-body. The most important conclusion from 
this figure is that detection of the charged Higgs bosons will be possible at future e + e~ 
colliders. 

E. Total Production of the Higgs Bosons 

We present in Tables I, II our results for the total production of h°, H°, A , H ± Higgs 
bosons, taking different values of the center-of-mass energy y/s, fundamental supersymmetry 
parameter tan (3, luminosity £, and the mass of the pseudoscalar m^o. We take the following 
values representative of the supersymmetric parameters m A o = 100 GeV and tan/3 = 10, 30. 
From these results we observed a strong dependence of the supersymmetric parameter tan /?, 
in particularly for tan (5 large, as well as, of the center-of-mass energy y/s, in the production of 
the different Higgs bosons (h°, H°, A , i? ± ). These results make apparent the big importance 
of investigate the possibility of detecting the Higgs bosons with the reactions at three- 
body e + e~ — > bbh°(H°,A ) and e + e~ — > t~v t H + ,t + v t H~ , at next generation linear e + e~ 
colliders. 
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Total Production of Higgs 


Bosons 


C= 


=50 fb' 1 


tan/3 = 


10 


tan/3 = 30 


h° 




1600 




1800 


H° 




700 




470 


A 




900 




935 


H+H- 




7000 




6500 



Table I. Total production of Higgs bosons h°, H°, A , F± of the MSSM for yfe = 500 GeV 

and m A o = 100 GeV. 



Total Production of Hig 


gs Bosons 


£=100 fb' 1 


£=200 fb' 1 


tan/3 = 10 


tan /3 = 30 


tan/3 = 10 


tan (3 = 30 


h° 




960 


1150 


1920 


2300 


H° 




370 


220 


740 


440 


A 




560 


615 


1120 


1230 


H+H- 




4700 


4900 


9400 


9800 



Table II. Total production of Higgs bosons h°, H°, A , #± of the MSSM for = 1 TeV 

and m A a = 100 GeV . 

IV. CONCLUSIONS 

In this paper, we have calculated the production of the neutral and charged Higgs bosons 
in association with 6-quarks and with ru T leptons via the processes e + e~ — > bbh, h = 
h°,H°,A° and e + e~ — > t~v t H + ,t + v t H- using the helicity formalism. We find that these 
processes could help to detect the possible neutral and charged Higgs bosons at energies of 
a possible Next Linear e + e~ Collider when tan/3 is large. 

In summary, we conclude that the possibilities of detecting or excluding the neutral and 
charged Higgs bosons (h°, H°, A , H^) of the Minimal Supersymmetric Standard Model in 
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the processes e + e~ — > bbh, h = h°,H°,A° and e + e~ — > t~v t H + ,t + is t H~ are important 
and in some cases are compared favorably with the dominant mode e + e~ — > (A , Z°) + /i, 
ft, = ft , A and e + e" — > H + H~ with tan/5 large. The detection of the Higgs boson will 
require the use of a future high energy machine like the Next Linear e + e~ Collider. 
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FIGURE CAPTIONS 



Fig. 1 Feynman Diagrams at tree-level for e + e — > bbh . For e + e — > one has to 

make only the change sin a/ cos/5 — > cos a/ cos/3. 

Fig. 2 Feynman Diagrams at tree-level for e + e~ — > 66 A . 

Fig. 3 Feynman Diagrams at tree-level for e + e _ — > t~v t H + ,t + u t H~ . 

Fig. 4 Total Higgs production cross sections e + e~ — > and e + e~ — > (A , Z°) + /i° with 
v^s = 500 and 4 different values of tan/3, 6, 10, 30, 50. We have taken m t = 175 GeV and 
M~ = 500 GeV and neglected squark mixing. 

Fig. 5 Same as in Fig. 4, but for e + e~ -> teii and e + e~ -> (A , Z°) + if . 

Fig. 6 Same as in Fig. 4, but for e+e" -> bbA° and e+e" -> H°) + A . 

Fig. 7 Same as in Fig. 4, but for e + e~ — > t~v t H + , t + v t H~ and e + e~ — > H + H . 
Fig. 8 Total Higgs production cross sections for ^/s = 1 T"eV and 4 different values of tan /3, 
6, 10, 30, 50. We have taken m t = 175 GeV, M~ = 500 GeV and neglected squark mixing. 
We display contours for e + e~ — > and e + e~ — > (A , Z°) + in the parameters space 
(<7 T - m A o). 

Fig. 9 Same as in Fig. 8, but for e + e~ -> M° and e+e" -> (A , Z°) + if . 
Fig. 10 Same as in Fig. 8, but for e+e~ -> 66A° and e+e~ -> if ) + A . 
Fig. 11 Same as in Fig. 8, but for e + e~ — > t~u t H + ,t + v t H~ and e + e~ — > H + H~. 
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